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ABSTRACT 


The  modern  maritime  battlefield  is  dominated  by  the  new 
generation  of  sea-skimming,  high-speed,  stealthy  and  highly 
agile  anti-ship  missiles.  Anti-ship  cruise  missile 
technology  continues  to  evolve  and  to  overcome  the 
performance  of  the  existing  ship  self-defense  weapon 
systems.  The  Free  Electron  Laser  (FEL)  could  be  the 
ultimate  speed-of-light ,  hard-kill  weapon  system,  offering 
unique  features  such  as  tunability,  high  power,  pinpoint 
accuracy  and  infinite  magazine.  Multimode  computer 
simulations  were  used  to  explore  the  operation  of  the  Thomas 
Jefferson  National  Acceleration  Facility  (TJNAF)  FEL  with 
untapered  and  positively  tapered  undulator.  The  final 
steady-state  power,  the  steady-state  gain  and  the  electron 
energy  spread  as  a  function  of  desynchronism  were  determined 
for  both  34.5  Mev  and  47.5  Mev  electron  beam  energies. 

This  thesis  also  includes  an  experimental  study  of 
damage  induced  to  Polyimide  Fiberglass  and  F2  Epoxy  samples, 
by  the  TJNAF  FEL.  Irradiations  of  the  samples  were 
conducted  changing  various  parameters,  such  as  the 
wavelength,  average  power,  pulse  repetition  frequency,  cross 
wind  and  spot  size  in  order  to  explore  the  damage  mechanism. 
At  this  stage  of  evolution,  TJNAF  FEL  is  now  capable  of  500W 
output  average  power.  In  order  to  achieve  the  required 
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to  a  small 


intensity  of  10  kW  /  cm2  the  beam  was  focused 
radius,  and  scaling  guidelines  were  developed 
predict  the  damage  caused  by  a  high  power  laser 
area . 


in  order  to 
over  a  large 
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EXECUTIVE  SUMMARY 


The  Free  Electron  Laser  (FEL)  could  be  the  solution  to 
the  complex  problem  of  ship's  self  defense  against  modern 
anti-ship  missiles  (ASM).  The  FEL's  immediate  response, 
infinite  magazine,  speed-of-light  target  engagement,  and 
continuous  tunability  make  it  an  excellent  alternative  of 
the  current  Close  In  Weapon  Systems  (CIWS) . 

Computer  simulations  were  used  to  explore  the  operation 
of  the  Thomas  Jefferson  National  Acceleration  Facility 
(TJNAF  FEL)  with  untapered  and  positively  tapered 
undulators.  Taper  rates  of  5=+4n,  +6n,  and  +8n  were  used. 
The  final  steady-state  power,  the  steady-state  gain  in  weak 
fields,  and  the  electron  energy  spread  as  a  function  of 
desynchronism  were  determined  for  both  34.5  Mev  and  47.5  MeV 
electron  beam  energy.  The  shape  of  those  curves  provides 
useful  information  about  the  FEL  operation.  The  final 
power  and  gain  were  significantly  smaller  for  the  47.5  MeV 
electron  beam  than  for  the  34.5  MeV  electron  beam,  because 
of  the  lower  value  of  dimensionless  current  j .  However,  the 
general  shape  of  the  curves  is  very  similar  for  each  energy. 
The  desynchronism  curves  present  a  sharp  peak  at  small 
values  of  desynchronism  accompanied  by  the  trapped-particle 
instability.  The  gain  versus  desynchronism  curves  peak  at 
moderate  values  of  d  and  give  near  few  gain  for  both  small 


xi 


and  large  values  of  desynchronism.  Limit  cycle  behavior  was 
observed  in  the  final  power  for  both  the  34.5  MeV  and  47.5 
MeV  electron  beams,  but  not  for  all  tapers.  The  steady- 
state  power  and  the  steady-state  gain  turned  out  to  be 
higher  for  the  conventional  undulator  than  for  the  tapered 
undulator.  Furthermore,  the  tapered  undulator  induces  more 
energy  spread  to  the  electron  beam,  and  demonstrates  less 
efficiency  for  beam  recirculation.  The  width  of  the 
desynchronism  curve  decreased  as  the  tapering  rate  is 
increased.  The  operating  width  of  the  desynchronism  curves. 
Ad,  as  a  function  of  taper  has  been  measured  in  experiments 
conducted  with  the  IR  Demo  FEL  at  Jefferson  Lab  [Ref  16]  . 
Simulation  results  are  summarized  in  Table  1.  There  is  good 
agreement  between  theory  and  experiment. 


Electron  Beam 

Energy 

34.5  MeV 

47.5  MeV 

Taper-  5 

0 

6n 

8n 

0 

4n 

6n 

8n 

Maximum  Power 

(dimensionless) 

206 

170 

112 

89 

9 

Maximum 

Gain  ( % ) 

85 

58 

33 

12 

76 

53 

32 

12 

De synchronism 
Curve 

Width  ( 

x  NA) 

0.38 

0.3 

0.19 

0.06 

0.34 

0.26 

0.18 

0.04 

Limit  Cycle 
Behavior 

YES 

NO 

NO 

YES 

YES 

NO 

YES 

Maximum 

Electron  Energy 
Spread  ( % ) 

8 

8 

7.5 

7 

7 

6.5 

6 

2 

Table  1.  Summarized  Simulation  Results. 
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This  thesis  also  includes  an  experimental  study  of 
damage  induced  to  missile  nose-cone  materials  by  the  FEL. 
Experiments  were  conducted  with  the  TJNAF  FEL  in  order  to 
study  the  damage  caused  by  the  laser  beam  to  some  sample 
materials  commonly  used  in  missile  cones.  The  TJNAF  FEL  is 
capable  at  this  point  of  500  W  output  power.  Since  1  MW 
power  over  a  spot  of  100  cm2  is  approximately  required  to 
shoot  down  a  missile,  the  beam  was  focused  to  a  small  spot 
area  in  order  to  achieve  the  required  intensity  of  10 

kW  /  cm  .  In  this  way  several  guidelines  were  developed  that 
could  predict  the  damage  caused  by  a  high  power  laser  over  a 
larger  area.  The  pulse  repetition  frequency  (PRF) ,  the 
wavelength,  the  average  power,  the  spot  size  and  the 
crosswind  speed  were  varied  during  the  irradiations  in  order 
to  determine  the  influence  of  these  parameters  on  the  damage 
induced.  Airflow  of  60-80  mph  speed  did  not  have  a 
significant  effect,  and  increased  the  penetration  depth  rate 
by  10%-15%.  The  penetration  depth  rate  is  proportional  to 
the  laser  intensity.  The  A=3.1  pm  wavelength  was  more 
effective  than  the  A=4.875  pm  wavelength.  In  addition,  the 
lower  PRF  (18.7  MHz)  seems  to  contribute  to  higher 
penetration  rates  due  to  the  higher  energy  per  pulse. 
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I. 


INTRODUCTION 


In  recent  years  highly  sophisticated  anti-ship  missiles 
(ASMs)  have  been  introduced  to  the  sea  warfare  environment 
making  the  modern  maritime  battlefield  much  more  compressed 
and  deadlier  than  ever  before.  The  ability  to  react  quickly 
and  correctly  to  multiple  ASM  attacks  is  now  an  essential 
part  of  the  self-defense  capability  of  a  warship.  Current 
ship  weapon  systems  tend  to  lose  their  ability  to  meet  this 
requirement  as  the  offensive  missile  technology  evolves 
spectacularly.  A  Free  Electron  Laser  (FEL)  ,  as  a  speed  of 
light,  hard-kill  weapon  system,  is  a  potential  solution 
offering  many  attractive  characteristics  in  performance, 
readiness,  and  life  cycle  management. 

Chapter  II  of  this  thesis  describes  the  FEL,  and  its 
principles  of  operation,  and  gives  an  overview  of  background 
theory.  Chapter  III  addresses  the  benefits  of  the  FEL  as  a 
hard  kill  weapon  system  against  anti-ship  missiles,  and 
analyzes  the  shortcomings  and  needs  of  current  ship  self 
defense. 

Chapter  IV  describes  the  operation  of  the  tapered 
undulator  design,  and  simulates  the  performance  of  the 
Thomas  Jefferson  National  Acceleration  Facility  (TJNAF)  FEL. 
Multimode  simulations  are  used  to  analyze  the  operation  with 


1 


untapered  and  positively  tapered  undulator,  describing  the 
evolution  of  short  optical  pulses  at  far  infrared 
wavelengths.  Trapped  particle  instability  and  limit  cycle 
behavior  are  observed  and  analyzed.  Furthermore  the  steady 
state  power,  the  steady  state  gain  and  the  electron  energy 
spread  as  a  function  of  desynchronism  are  determined  for 
both  34.5  Mev  and  47.5  Mev  electron  beam  energies. 

Chapter  V  presents  the  analysis  of  experimental 
irradiations  on  several  materials  commonly  used  in  missile 
bodies.  The  goal  of  the  experiments  was  to  determine  the 
amount  of  damage  induced  by  the  unique  FEL  beam  of  short 
pulses  and  develop  scaling  laws  that  could  predict  the 
damage  caused  by  a  much  higher  laser  power  over  a  large 
area . 


2 


II. 


FREE  ELECTRON  LASER  THEORY 


A.  DESCRIPTION  AND  BASIC  OPERATION 

The  Free  Electron  Laser  (FEL)  concept  was  first 
introduced  by  John  Madey  in  1971.  The  basic  parts  of  a  FEL 
are  the  electron  accelerator,  the  undulator,  and  the  optical 
resonator  cavity  as  shown  in  Figure  1. 


The  undulator  is  a  configuration  of  magnets  or 
superconducting  coils,  which  produce  a  spatially  periodic 
magnetic  field.  A  relativistic  electron  beam  generated  by 
the  accelerator  is  injected  into  the  undulator.  As 
electrons  travel  through  the  undulator,  they  are 
periodically  deflected  in  the  transverse  direction  by  the 
imposed  magnetic  field.  This  transverse  motion  causes  the 
electrons  to  radiate  (spontaneous  emission)  in  a  narrow 
forward  cone.  The  electromagnetic  radiation  is  stored  in 
the  optical  cavity  and  the  optical  power  grows  over  many 
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passes  inside  the  resonator.  The  coupling  of  the  optical 
field,  the  magnetic  field,  and  the  electron  beam  allows 
extraction  of  kinetic  energy  from  the  relativistic  electrons 
and  conversion  to  high  optical  power.  Subsequent  electrons 
moving  inside  the  undulator  interact  with  the  continuously 
growing  optical  field  and  radiate  in  the  presence  of  light; 
this  is  stimulated  emission. 

The  combined  effect  of  the  undulator  magnetic  field  and 
the  electromagnetic  radiation  gives  rise  to  the  axial 
bunching  of  the  electrons  at  the  optical  wavelength, 
allowing  them  to  radiate  coherently. 


B.  PARAMETERS  AND  TYPICAL  VALUES 

The  FEL's  most  important  operation  parameters  are 
listed  in  Table  2. 


PARAMETER 

DEFINITION 

RANGE 

B 

Undulator  Magnetic  Field  (rms) 

2-7  KG 

N 

Number  of  Undulator  periods 

~  100 

Ao 

Undulator  Wavelength 

2-10  cm 

A 

Optical  Wavelength 

5nm-5cm 

L 

Undulator  Length 

1-20  m 

K 

Undulator  Parameter 

0.1-10 

Y 

Lorentz  Factor 

5-5000 

vmc2 

Electron  Beam  Energy 

MeV-GeV 

I 

Electron  Beam  Current 

1-1000A 

rb 

Electron  Beam  Radius 

!  1mm- 1cm 

Pin 

Input  Power 

10-100  MW 

P  out 

Output  Power 

1-10  MW 

n 

Efficiency  of  Energy  transfer 

l%-50% 

Table  2.  FEL's  Operation  Parameters. 
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The  undulator  parameter  K  is  defined  by  K  =  eBA0  /  2nmc2 
and  its  value  is  usually  around  1.  The  radiation  optical 
wavelength  is  A  =  A0(l  +  K2)/2y2  (Ref.  1).  The  FEL's 

operating  wavelength  can  be  easily  tuned  simply  by  adjusting 
the  electron  beam  energy  or  the  undulator  parameter.  It  is 
worth  noting  that  the  FEL  can  be  designed  to  work  in  an 
extremely  wide  spectrum  of  wavelengths  and  in  new  wavelength 
ranges,  where  there  are  few  powerful  sources  of  radiation 
like  20pm  to  500pm  and  20nm  to  300nm.  This  remarkable 
feature  makes  the  FEL  an  attractive  device,  demonstrating 
its  potential  for  various  applications. 


C.  ELECTRON  EVOLUTION  INSIDE  THE  UNDULATOR 

1 .  The  Undulator 

Two  kinds  of  undulator  polarizations  are  commonly  used 
in  FELs,  the  helical  and  the  planar.  The  helical  undulator 
consists  of  superconducting  coils  and  generates  a  magnetic 
field  B  =  B  [cos(k0z)  ,  sin(k0z)  ,  0]  on  axis  with  k0  =  2n  /  A0 , 

and  it  produces  circularly  polarized  light.  The  planar 
undulator  consists  of  an  arrangement  of  magnets  and 
generates  a  magnetic  field  B  =  B  [0  ,  sin(k0z)  ,  0]  on  axis, 
producing  linearly  polarized  light. 
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Figure  2.  Helical  Undulator.  Figure  3.  Linear  Undulator. 

(From  Ref.  2)  .  (From  Ref.  2)  . 


2 .  Resonance  Condition 

As  electrons  move  through  the  undulator  in  the  axial 
direction,  they  are  periodically  deflected  in  the  transverse 
direction  by  the  spatially  varying  magnetic  field.  A 
typical  motion  of  electrons  inside  the  undulator  is  shown  in 
the  figure  below.  Note  that  z  and  x  are  the  axial  and 
transverse  directions  respectively. 


Figure  4.  Electron  motion  inside  the  Undulator. 
(From  Ref.  2) . 
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As  an  electron  moves  through  an  undulator  wavelength 
X0,  it  completes  a  full  oscillation  period  in  the  transverse 
direction  and  emits  a  wavelength  of  light  X.  So  one 
wavelength  of  light  traveling  with  speed  c  passes  over  a 
slower  electron,  as  it  travels  through  an  undulator 
wavelength  A0,  with  speed  (3zc .  This  particular  electron- 
photon  "race"  is  called  the  resonance  condition  and  is  shown 
schematically  in  the  figure  below. 


Figure  5.  Electron  -  Photon  Race.  (After  Ref.  3). 


The  electron  acceleration  causes  the  emission  of 
synchrotron  radiation  in  a  cone  of  width  1/y  (Ref.  2), 

where  y  =  1/^1  -  (B2  is  the  Lorentz  Factor  (3  =  u/c,  and  u  is 


the  electron's  speed. 


We  define  (3zc  the  electron  speed  in  the  axial 


direction,  where  c  is  the  speed  of  light.  We  can  easily 
calculate  the  wavelength  A  of  the  emitted  radiation  using 
the  electron-photon  speed  difference. 

f  \ 


A  =  (c  -  p2c) 


(\  \ 

=  ^  •  (i  -  a, 

:)  =  A0 

M 

3z 

^0 


2y2  -  1 


3- 


For  highly  relativistic  electrons,  y  »  1  and 

A  *  A0  /  2y2  . 

3 .  Pendulum  Equation 

We  consider  a  helical  undulator  magnetic  field 
B  =  B  [cos(k0z)  ,  sin(k0z)  ,  o],  where  k0  =  2n/A0  is  the 

undulator  wave  number  and  A0  is  the  undulator  wavelength. 
We  assume  for  the  optical  field  a  circularly  polarized  plane 
wave,  with  electric  and  magnetic  fields 

Es  =  E  [cosT  ,  -sinT  ,  o]  and  Bs  =  E  [sinT  ,  cosT  ,  o]  respectively 
where  T  =  kz-«t  +  cp,  (p  is  the  optical  phase, 

k  =  2n/A  is  the  optical  wave  number,  and  A  is  the  optical 
wavelength.  It  should  be  pointed  out  that  the  electric  and 
magnetic  fields  Es,  Bs  have  the  same  amplitude  E,  as  they  are 
measured  in  cgs  units. 
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Relativistic  electrons  with  charge  magnitude  e  are 
injected  in  the  undulator  with  energy  ymc2  and  momentum 
P  =  ymu ,  where  y  is  the  Lorentz  factor,  m  is  the  electron 

mass,  u  =  j3c  =  pxc  +  (3zc  is  the  electron's  velocity,  and 

P_i_c,  p2c  are  the  velocity  components  in  the  transverse  and 
axial  direction  respectively.  Electrons  interact  with  the 
fields  and  they  experience  relativistic  force 


dP 

F  =  —  =  -e| 
dt 


-u- 
E  +  —  x  Bt 


[1] 


From  [1] 


d(yP)  _  ^_e  fe 

dt  me 


[E  +  P  x  Bt], 


[2] 


where  BT  is  the  total  magnetic  field  of  the  undulator  and 
optical  field  combined  together  BT  =  B  +  Bs .  The  magnetic 
field  is  always  perpendicular  to  the  electron  motion  and 
does  not  affect  the  electron's  energy  ymc2  .  The  electron 
energy  change  is  given  by 


d(ymc  )  -  -  _  dy  e  -r  - 

— - -  =  F  •  u  =  -eE  •  u  =>  —  = - p  •  E  . 

dt  dt  me 

Inserting  the  fields  into  the  Lorentz  equation  we  get 


[3] 


=  — -  [e  ( 1-  3Z)  •  (cos  -  sin  T,  0  )  +  32B  •  (-  sin  ( koZ ) ,  cos  ( k0  z ) ,  O)]  [  4  ] 

dt  me 

—  ^  =  — -  [o,  0,  e((3xcosW  -  PySin1?)  +  B  ■  (pxsin(k0z)  -  (3y  cos(k0z))]  [5] 
dt  me 
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~z  =  -  —  e[Px  cos  ^  ~  Py  sin  w].  [6] 

dt  me  * 

For  relativistic  electrons  E(1  -  pz)  «  BPZ,  so  we  can 
disregard  the  term  E  { 1  -  pz)  •  (cos1?  ,  -  sin1?  ,  0) .  Within  this 
approximation,  integration  of  equation  [4]  with  respect  to 
time  gives  the  transverse  motion  as 


Pi 


K  r 

- [cos  (k0z) 

Y 


sin  (kgz)  ,  0] 


[7] 


where  K  is  the  undulator  parameter  and 


Px 

Py 


K 

- cos  (k0z) , 

Y 


- sin  (k0z) . 

Y 


[8] 

[9] 


We  substitute  Px, Py  into  equation  [6]  and  obtain 


dy  .  eKE  ,  v 

—  =  Y  =  - cos(£  +  cp). 

dt  ymc 


[10] 


The  Lorentz  Factor  y  describes  the  electron's  energy, 
£  =  (k  +  k0)  z  -  cot  is  the  electron's  phase  generated  by  the 
combination  of  the  undulator  and  optical  fields,  and  <j>  is 
the  optical  phase.  Since  A  «  A0  ,  k  »  k0  and  y»l 

(relativistic  electrons),  the  initial  electron  phase  in  the 
beginning  of  the  undulator  at  t=0  is  =  (k  +  k0)z0  «  kz0 . 
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Combining  this  with 


From  equation  [7],  we  have  p2 


K 


the  Lorentz  factor,  we  get 


y“2  _  2  _  ^  =  ]__P2_  =  ]_  —  P2  —  —  «  1  —  P2 

y2 


[ii] 


Differentiation  of  [11]  with  respect  to  time  gives 


1  =  1  M* 

Y  (l  +  K2) 


[12] 


which  relates  y  to  P2 .  Differentiating  £  twice,  we  acquire 
t,  =  (k  +  k0)z  -  to  =  (k  +  k0)cPz  -  co  [13] 


and 

?  =  (k  +  k0)cp  =>  pz  =  - — - . 

[14] 

(k  +  k0)c 

We 

substitute  equation  [14]  into 

[12] 

and  get 

i  = 

7i - n - • 

[15] 

Y 

(l  +  K2j-  (k  +  k0)c 

We  now  make  use  of  the  resonance 

condition 

X  = 

X0(l  +  K2)  /  2y2  or  k0  =  k(l  +  K2)  /  2y2 ,  and 

assume 

that  for 

relativistic  electrons  Pz  =  1  and  k  +  k0  = 

k,  since  k  »  k0  . 

Equation  [15]  becomes 

i  = 

l 

[16] 

Y 

2k0c 

Combining  equations  [12]  and  [16],  and  solving  for  C,  we 
obtain 
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[17] 


W  2k0eKE  ,  v 

C  =  - 2 -  COSte  +  9)  , 

y  m 

which  describes  the  dynamics  of  the  electron  phase  by  the 

Pendulum  Equation. 

The  time  interval  for  an  average  electron  to  cross  the 
undulator  is  At  =  L/(30c,  which  motivates  the  introduction 

of  the  relevant  dimensionless  time  i  =  (30ct  /  L  =  ct  /  L  , 
ranging  from  0  to  1  along  the  whole  undulator  length  L.  We 


00 

use  the  notation  of  £  as  the  second  derivative  of  £  with 

respect  to  x,  and  we  define  |a|  =  4nNeKLE  /  y2mc2  as  the 

dimensionless  optical  field  amplitude.  The  pendulum 

equation  takes  its  final  form 
00 

£  =  M  cos  (  £  +  cp).  [18] 

Certain  important  remarks  must  be  made  from  this  analysis: 

•  C,  is  a  microscopic  variable  describing  the 
relative  position  of  an  electron  within  an 
optical  wavelength,  as  it  travels  through  the 
undulator. 

•  Equation  [10]  relates  y  ,  which  is  the  electron's 
rate  of  energy  change,  to  the  electron's  phase 
(£  +  <p).  As  cos(£  +  cp)  takes  values  from  -1  to  1 
within  one  optical  wavelength,  half  of  electron's 


12 


phases  result  in  y  >  0,  and  the  other  half  give 
y  <  0,  which  implies  that  electrons  become  spread 
in  energy. 

•  The  fact  that  y  °c  EK  /  y  implies  that  stronger 
fields  result  in  a  larger  electron  energy  change, 
whereas  more  energetic  electrons  (larger  y) 

reduce  the  coupling  of  the  undulator  and  optical 
fields . 

•  Electrons  becoming  spread  in  energy  within  an 
optical  wavelength,  implies  that  some  of  them 
gain  energy  and  move  faster  while  other  lose 
energy  and  move  slower,  thus  they  become  bunched 
on  the  scale  of  the  optical  wavelength. 

•  Electron  motion  inside  the  undulator,  which  is 
rather  complicated,  is  described  by  the  well- 
known  pendulum  equation ,  which  facilitates 
remarkably  the  visualization  of  the  electron 
phase  dynamics,  and  the  understanding  of  the 
phase  space  diagrams. 
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D .  DIMENSIONLESS  PARAMETERS 

Several  dimensionless  variables  can  summarize  recurring 
combinations  of  physical  parameters,  which  govern  the  FEL 
operation,  and  are  used  to  gain  insight  into  the  relevant 
physical  processes.  An  important  quantity  characterizing 
many  FEL  properties  is  the  undulator  parameter 

K  =  eBA0  /  2nmc2  where  B  is  the  rms  field  strength  over  each 
period,  A0  is  the  undulator  wavelength,  and  e,  m,  are  the 

electron  charge  and  mass  respectively.  In  most  FELs  K  «»1. 
A  pass  of  the  electron  pulse  through  the  undulator  is 
described  by  the  dimensionless  time  x  =  ct  /  L  where 
t  =  0  — »  1  .  The  electron  phase  £  and  phase  velocity 

v  =  ^  =  L[(k  +  k0)  •  -  k],  describe  the  electron  velocity  and 

position  in  a  section  of  the  beam  one  wavelength  long,  and 
define  bunching  and  bunching  rate  of  the  electrons.  The 
phase  and  phase  velocity  are  coupled  by  the  dimensionless 

complex  optical  field  a  =  |a|e1<p  through  the  simple  pendulum 

OO 

equation  £  =  |a|  cos  (  £  +  <p) .  The  field  amplitude  is 

|ot|  =  4nNeKLE  /  y2mc2 ,  and  determines  the  electron-bunching 
rate  on  the  optical  wavelength  scale.  When  a  «  n  the 
optical  field  is  weak,  and  when  a  »  n  the  FEL  operates  in 
the  strong  field  regime.  The  response  of 
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the  optical  wave  to  bunching  in  the  beam  is  described 
by  the  dimensionless  current  j  =  8N(enKL)2ne  /  y3mc2 ,  where  ne 
is  the  electron  density  within  the  electron  beam.  The 
dimensionless  current  provides  the  coupling  between  the 
electron  beam  and  the  light  wave  through  the  wave  equation 

a  =  -j(e  y  (Ref.  1)  .  When  j«n  the  FEL  gain  is  low,  and 
when  j»n  the  FEL  gain  is  high.  Each  definition  above,  £ 
and  v  describing  the  electrons,  a  describing  the  light 
wave,  and  j  describing  the  whole  FEL,  is  physically 
meaningful  and  useful  in  evaluation  of  many  FEL  effects. 

E .  GAIN  SPECTUM 

Figure  6  shows  the  FEL  gain  spectrum  for  a  typical  low 
current  FEL  with  j=2  and  N=21,  as  the  initial  optical  field 
amplitude  a0  evolves  from  0  to  40.  G  is  the  FEL  gain  for 
one  pass  of  the  electron  beam  through  the  undulator  in  a 
single  optical  mode,  and  v0  is  the  initial  electron  phase 
velocity  at  the  beginning  of  the  undulator.  Figures  7,  8 

and  9  present  a  cross  section  of  the  gain  spectrum  for 
a0  =  0,  a0  =  20  and  a0  =  40  respectively.  The  plots  were 
generated  by  computer  simulations.  The  gain  has  a  maximum 
value  (peak  gain),  for  each  optical  field  strength  a0 . 

Furthermore,  the  position  of  the  peak  gain  changes  as  a0 
increases,  and  the  gain  spectrum  distorts. 
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Figure  6.  FEL  Gain  Spectrum,  j=2. 


Figure  7.  FEL  Gain  Spectrum,  j=2,  a0 


0  . 


Figure  8.  FEL  Gain  Spectrum,  j=2,  a0  =  20. 


The  value  and  the  position  of  the  peak  gain  are  crucial 
to  the  FEL  operation,  since  they  determine  the  optical 
frequencies  that  are  going  to  be  amplified  as  well  as  the 
final  optical  power.  The  FEL  gain  is  clearly  antisymetric 

about  resonance  vo=0.  In  weak  optical  fields  the  gain 
spectrum  has  a  peak  value  of  G  «  0.13j  at  vo  =  2 . 6 ,  and 
width  Av0«n.  The  gain  is  zero  at  resonance.  As  the 
optical  field  grows  to  saturation  the  gain  spectrum 
distorts,  and  becomes  broader  with  width  Av0  « •  The  peak 
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gain  decreases  substantially  and  moves  to  higher  values  of 
v0 .  In  the  case  of  j  =  2,  the  peak  gain  drops  from  27%  in 
weak  optical  fields,  to  1%  when  the  optical  field  saturates. 
The  phase  velocity  of  peak  gain  moves  from  vc=2.6  to 

vo=10.8.  The  corresponding  change  in  the  wavelength  is 
AX/A  =  Av0/4nN  =  3%  in  this  example. 

F.  FEL'S  APPLICATIONS 

The  FEL  mechanism  can  be  designed  to  operate  over  a 
large  region  of  wavelengths,  from  nanometers  to  centimeters. 
Each  FEL  also  has  the  capability  of  being  continuously  tuned 
simply  by  changing  the  electron  beam  energy.  Most  often  the 
FEL  is  operated  so  that  the  net  energy  transfer  is  from  the 
electron  beam  to  the  radiation  field  in  order  to  amplify  it. 
Nevertheless,  it  is  possible  to  operate  the  FEL  so  that  the 
net  energy  transfer  is  from  the  radiation  field  to  the 
electrons  making  an  accelerator  or  an  inverse  FEL. 

The  FEL  has  become  an  exciting  alternative  to  other 
radiation  sources,  like  microwave  tubes  and  lasers,  and  can 
extend  the  operational  range  of  both.  Because  of  its 
flexibility,  the  FEL  can  be  applied  in  many  diverse  areas 
such  as  lithography,  plasma  heating,  particle  acceleration, 
as  well  as  material,  biological,  medical,  and  solid  state 
research. 
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III.  THE  FEL  AS  SHIP  SELF-  DEFENSE  WEAPON 


A.  THE  ANTI -SHIP  MISSILE  THREAT 

In  the  modern  naval  battlefield,  anti-ship  missiles 
(ASM)  have  become  a  very  serious  threat  even  to  the  most 
modern  and  complex  warships.  Faster  and  faster  ASMs  have 
dramatically  decreased  response  times  for  surface  vessels. 
A  simple  scenario  helps  us  realize  the  complexity  of  the 
problem. 

Consider  a  sea  skimming  ASM,  flying  less  than  10  meters 
above  sea  level  and  heading  towards  a  ship.  Let  the  height 
of  the  ship's  sensors  be  25  meters  above  the  sea  surface,  so 
that  the  horizon  distance  in  which  the  missile  can  be 
detected  and  tracked  is  approximately  10  nautical  miles 
(nm) .  Given  the  missile's  speed  to  be  1  mach  =  600  knots  = 
10  nm/min,  it  takes  one  minute  for  the  missile  to  impact  the 
ship.  The  time  is  decreased  to  30  seconds  if  the  missile 
approaches  with  speed  mach  2. 

The  Aegean  Sea,  where  the  Hellenic  naval  combatants 
usually  operate,  has  a  very  specific  geography.  It  is  a 
complex  of  thousands  of  neighboring  islands,  less  than  12 
nautical  miles  apart  each  other.  This  complicated  battle 
arena  gives  rise  to  certain  specific  scenarios.  For 
instance,  a  battleship  operating  5  nautical  miles  from  the 
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coast  of  a  small  island  could  suddenly  be  face  to  face  with 
a  Harpoon  ASM  emerging  behind  the  island,  using  the  waypoint 
navigation  feature. 

Due  to  the  spectacular  development  of  anti-ship  missile 
technology,  ASMs  have  become  state-of-the  art  weapons  with 
highly  supersonic  speeds,  ten's  of  g's  maneuverability, 
multiple  guidance  modes,  countermeasure  capability,  reduced 
radar  cross  section  and  IR  signature,  and  pinpoint  accuracy. 
There  seem  to  be  no  end  to  this  development  making  it  hard 
for  the  ship  self-defense  systems  to  keep  pace. 


B.  SHIP  SELF-DEFENSE 

The  existing  shipboard  anti-missile  systems  are  divided 
into  two  major  categories:  passive  and  active.  Passive 
methods  include  decoys,  chaffs,  stealth  technology,  as  well 
as  tactics  and  techniques,  which  alter  the  ship's  IR  or 
RADAR  signature,  and  generally  deceive  and  confuse  the 
missile  sensors. 

On  the  other  hand,  active  systems  imply  strong 
interaction  with  the  incoming  missile  target  and  they  are 
further  divided  into  two  categories:  soft  kill  and  hard 
kill.  Soft  kill  methods  disable  the  threat  missile  by 
attacking  its  electronics  or  blinding  its  sensor,  but 
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otherwise  leaving  the  target  body  undamaged.  Hard  kill 
weapon  systems  actually  destroy  the  target-missile  achieving 
large-scale  physical  damage.  We  are  going  to  be  focused  on 
the  hard  kill  method,  which  is  the  most  desirable.  A 
scenario  of  layered  ship  self-defense,  involving  both  soft 
kill  and  hard  kill  methods,  is  demonstrated  in  Figure  10. 
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C.  AVAILABLE  HARD  KILL  WEAPON  SYSTEMS 

Several  hard  kill  weapon  systems  are  used  to  confront 
the  complicated  ASM  threat. 

1 .  Large  Caliber  Guns 

The  5/54  or  3/62  OTOMELARA  guns  have  effective  ranges 
of  10  and  5  nautical  miles  respectively,  and  are  the 
standard  Hellenic  Navy  guns  that  can  be  used  as  a  medium 
range  defense,  but  with  extremely  low  efficiency.  The  fire 
control  system,  the  shell's  speed,  and  the  relatively  slow 
rate  of  fire  are  not  capable  of  shooting  down  an  approaching 
supersonic  sea  skimming  ASM  with  high  maneuverability. 
Moreover,  the  shell  dispersion  (increasing  with  large  range) 
limits  the  accuracy  considerably  and  lowers  the  probability 
kill  to  minimal  values. 

2 .  Long  Range  Surface-Air  Missiles 

Another  approach  considered  as  a  long-range  defense  is 
the  use  of  antimissile  missiles.  Several  surface-to-air 
missiles,  such  as  the  standard  Navy  missile  (SM-1,  SM-2)  or 
the  Nato  Sea  Sparrow  (NSSM)  have  been  developed,  and  are 
equipped  with  excellent  flight  characteristics, 
maneuverability,  and  effective  range. 

However,  in  order  to  take  advantage  of  the  long 
effective  range,  a  long  detection  range  is  also  required. 
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Since  the  sea  skimming  ASM  missile's  flight  altitude  is 
about  8  m  with  a  corresponding  detection  range  at  the 
horizon  of  11  nautical  miles,  making  the  long-range  missile 
advantage  becomes  useless.  Moreover,  physical  dynamics 
require  that  defensive  missiles  should  have  2  to  3  times  the 
maneuverability  of  the  incoming  missile.  Current  anti-ship 
missiles  have  the  ability  to  fly  at  supersonic  speeds, 
pulling  more  than  lOg  in  terminal  maneuvers.  A  simple 
calculation  makes  the  maneuvering  requirements  of  the 
defensive  missile  unfeasible. 

3.  CIWS 

The  innermost  layer  of  defense  against  anti-ship 
missiles  is  provided  by  the  Close  In  Weapon  Systems  (CIWS) . 
The  inner  layer  defense  systems  could  be  gun  based,  such  as 
Raytheon's  20mm  Phalanx,  FABA's  20mm  Meroka,  Signaal's  30mm 
Goalkeeper,  Oerlicon  -  Royal  Ordance  35mm  Millenium,  or 
OTOBreda ' s  twin  25mm  Myriad.  Some  Navies  advocate  short- 
range  missile  based  CIWS ' s  such  as  the  US-German  RAM  or  the 
French  Mistral.  The  basic  characteristics  and  features  of 
the  currently  most  used  CIWS  are  discussed  below. 
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The  Phalanx  Gun 


a; 

At  this  time,  the  Phalanx  Gun  is  the  standard  CIWS 
for  both  the  U.S  and  Hellenic  Navies.  It  is  the  most 
extensively  used  CIWS  in  the  world  with  more  than  700 
Phalanx  systems  deployed  by  21  Navies.  Since  1979,  when  it 
was  initially  deployed,  the  weapon  has  evolved  from  Block  0 
to  Block  IB  through  system  improvements,  upgrades,  and 
modifications.  The  system  is  a  rapid  fire,  computer 
controlled  radar  and  gun  system  that  automatically  carries 
out  search,  detection,  target  threat  evaluation,  tracking, 
and  firing.  It  is  equipped  with  an  electronically 
controlled,  pneumatically  driven,  20  mm  six-barrel  gun 
capable  of  firing  4500  rounds  per  minute  within  a  range  of  2 
km.  Although  the  weapon  provides  superior  performance 
against  anti-ship  missiles,  it  has  several  deficiencies.  It 
experiences  random  and  variable  bullet  dispersion,  which 
diminishes  its  ability  to  destroy  its  intended  targets  at  a 
range  where  missile  debris  does  not  hit  the  ship.  The  high 
rate  of  fire  causes  overheating  of  the  barrels,  so  that  the 
firing  must  be  limited  to  bursts  of  only  few  seconds.  The 
gun's  magazine,  which  is  capable  of  holding  up  to  1000 
rounds,  can  be  expended  in  just  four  to  five  engagements  and 
requires  up  to  30  minutes  to  be  reloaded.  If  the  missile 
attacks  are  more  frequent,  it  means  that  the  ship  is  left 
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with  no  self-defense.  Finally,  computer  simulations  have 
shown  that  the  probability  of  hitting  an  incoming  missile  is 
extremely  low  for  ranges  more  than  250  m  (Ref.  11)  .  In 
Figure  11  we  can  see  the  probability  of  a  single  phalanx  gun 
penetrator  hitting  an  incoming  missile  flying  at  altitude  of 
7  meters  above  sea  surface  with  speed  Mach  2,  versus  the 
missile  range  from  the  ship.  It  is  easy  to  see  that  the 
probability  drops  very  rapidly  as  the  range  increases. 


Figure  11.  Probability  of  a  single  round  to  hit  the  missile 
versus  missile  range. 
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Figure  12,  shows  how  the  number  of  hits  acquired 
by  a  missile  with  approach  speed  Mach  2  and  altitude  7  m 
above  the  sea  level,  varies  with  range  from  the  ship.  The 
dashed  line  is  for  a  firing  rate  of  3000  rounds/min  (Block 
1A  model)  ,  while  the  solid  line  is  for  4500  rounds/min 
(Block  IB  model)  .  Since  5  to  6  hits  are  needed  to  destroy 
the  missile  a  typical  killing  range  seems  to  be  around  150  m 
from  the  ship.  The  missile  may  be  destroyed  within  a  range 
of  100  to  200  m,  but  its  fragments  contain  enough  kinetic 
energy  to  reach  the  ship  and  cause  serious  damage. 


Figure  12.  Number  of  shots  acquired  by  the  missile  versus 
range . 
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Figure  13  presents  the  probability  of  a  missile 
fragment  hitting  the  ship  versus  the  missile  destruction 
range.  It  is  clear  that  the  probability  increases 
dramatically  when  the  missile  is  destroyed  closer  to  the 
ship.  For  a  missile,  approaching  with  Mach  2  speed  and  7  m 
altitude  destroyed  by  the  Phalanx  Gun  at  200  m  from  the  ship 
and  breaking  into  20  pieces,  there  will  be  about  10  of  them 
that  finally  hit  the  ship. 


Figure  13.  Probability  of  a  single  missile  fragment  to  hit 
the  ship  versus  the  missile  destruction  range. 
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b)  The  Rolling  Airframe  Missile  (RAM) 

Another  relatively  new  weapon  system  is  the 


American/German  Rolling  Airframe  missile  (RAM)  that  was 
first  introduced  in  1993  to  counter  the  anti-ship  cruise 
missile  threat.  It  is  a  lightweight  quick  reaction,  high 
firepower  guided  missile  weapon  system  that  provides  another 
layer  of  hard  kill  defense  and  fills  the  gap  between  the 
Phalanx  CIWS  and  the  longer-range  anti-missile  missile 
system  as  the  Nato  Sea  Sparrow. 

RAM  comes  with  a  21-round  guided  missile  launcher, 
which  provides  high  firepower.  The  supersonic  missile 
demonstrates  an  excellent  kill  performance  resulting  from 
its  high  maneuverability  and  accuracy  as  well  as  from  its 
effective  warhead  and  fusing.  Moreover,  the  latest  missile 
upgrade  RAM  Block  1  Infrared  (IRMV)  exhibits  multiple 
guidance  modes,  reduced  IR  and  Radar  signatures, 
countermeasures  capability,  and  improved  performance  in 
degraded  weather  conditions  such  as  heavy  rain  or  fog.  The 
currently  produced  RAM  missile  (Block  0)  has  been  fired  in 
105  flight  tests  resulting  in  100  successful  intercepts. 
Its  reliability  has  been  sustained  in  extensive  field  tests 
under  high  stress  conditions  and  realistic  scenarios  thus 
demonstrating  that  RAM  is  one  of  the  world's  most  advanced 
ship  self-defense  weapon  systems.  Another  system  that  is 
also  under  development  is  the  Raytheon's  SEA  RAM,  which 
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combines  elements  from  both  Phalanx  and  RAM  CIWS.  An  11- 


missile  round  RAM  launcher  assembly  loaded  with  RAM  B1 
guided  missiles  replaces  Phalanx's  20mm  gun.  SEA  RAM 
inherits  RAM's  accuracy,  extended  range  and  maneuverability 
combined  with  the  Phalanx's  Block  IB's  high  resolution 
target  detection  search  and  track  sensor  systems  and 
reliable  quick-response  capability.  SEA  RAM  is  an 
affordable  capability  upgrade  and  especially  attractive 
option  for  those  navies  that  have  already  deployed  the 
Phalanx,  like  the  Hellenic  Navy.  It  fits  the  exact 
footprint  of  the  Phalanx,  uses  the  same  power  and  requires  a 
minimal  ship  modification. 

Nevertheless,  anti-ship  missile  technology  evolves 
extremely  fast,  ASMs  have  not  yet  reached  the  peak  of  their 
evolution.  Future  high  performance  threats  are  expected  to 
be  stealthier,  faster,  and  lower  flying  missiles  able  to 
perform  unique  maneuvering  profiles  in  order  to  defeat  ship 
defenses.  Since  a  defensive  missile  is  required  to  have 
substantially  superior  speed,  guidance  system,  and 
maneuverability  than  the  incoming  threat,  it  is  a  matter  of 
time  that  RAM  or  any  current  missile  based  inner  layer 
defense  system  will  become  outdated  and  ineffective.  The 
defensive  -  offensive  missile  technology  race  is  clearly  won 
by  the  offensive  one. 
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D.  DEFENSE  AT  THE  SPEED  OF  LIGHT 

So  far,  only  kinetic  kill  methods  have  been  used  to 
project  power  to  the  target,  generally  in  the  form  of 
projectiles,  rockets  and  missiles.  Kinetic  energy  weapons 
require  a  finite  flight  time  to  intercept  a  target.  This 
time  turns  out  to  be  crucial  when  a  threat  is  approaching 
with  supersonic  speed,  making  the  reaction  time  for  the  ship 
minimum  and  leaving  no  chance  for  reengagement. 


Figure  14.  HEL  Engagement  Scenario.  (From  Ref.  7). 
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High  Energy  Lasers  (HELs)  are  now  being  developed  to  be 
used  as  weapons  delivering  lethal  power  at  the  speed  of 
light.  A  high  power  optical  beam  traveling  with  the  amazing 
speed  of  300,000  km/sec  is  now  available  to  boost  ship's 
self-defense.  The  HEL  provides  an  effective  engagement 
range  of  approximately  10  km  and  a  practically  zero  flight 
time  out  to  the  target.  In  the  case  of  a  typical  short- 
range  (10  km)  engagement  with  a  modern  defensive  missile 
system,  the  time  between  the  employment  and  threat 
destruction  can  be  as  big  as  10  seconds.  In  the 
unfortunate  case  of  a  miss  there  is  no  time  for 
reengagement.  On  the  other  hand  it  takes  only  30  jxsec  for 
the  laser  beam  to  reach  the  target  and  about  2  seconds  to 
destroy  it.  Such  a  short  time  provides  the  benefit  of  a 
"second  chance",  if  there  is  a  miss  and  moreover  the 
advantage  of  engaging  multiple  threats. 

An  example  of  an  HEL  is  MIRACL  (Mid- Infra-Red  Advanced 
Chemical  Laser) ,  which  has  been  developed  at  White  Sands 
Missile  Range  in  New  Mexico.  MIRACL  is  a  deuterium  fluoride 
laser  that  operates  at  3.8  pm  wavelength,  and  has  achieved 
an  output  power  of  2.2  MW.  At  tests  conducted  in  February 
of  1996  under  the  Nautilus  program,  MIRACL  successfully 
engaged  and  destroyed  a  short-range  rocket  in  flight  using 
the  Sea  Lite  pointing  and  tracking  system. 
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Undoubtedly,  the  HEL  could  enhance  substantially  a  ship's 
self-defense  capability,  and  HEL  development  should  be 
considered  seriously  to  provide  an  alternative  to  current 
short-range  defense  weapons.  It  is  the  weapon  of  the  21st 
century  and  will  always  be  effective.  Its  instantaneous 
response  and  amazing  capabilities  can  never  be  outdated,  as 
the  speed  of  light  will  never  be  exceeded  by  the  performance 
of  any  future  advanced  missile. 


E.  LASER  BEAM  LETHALITY 

The  HEL  weapon  concept  is  fundamentally  different  from 
that  of  the  kinetic  energy  weapon.  A  kinetic  energy  weapon 
requires  a  finite  time  to  deliver  its  warhead  or  mass  to  the 
target,  but  causes  damage  instantly.  On  the  other  hand,  the 
laser  beam  arrives  at  the  target  instantly,  but  it  takes  a 
finite  time  to  cause  damage.  The  damaging  effect  of  the  HEL 
is  a  thermal  process  nearly  identical  to  cutting  a  peace  of 
metal  with  a  blowtorch.  The  question  posed  is  how  lethal  a 
laser  beam  has  to  be  in  order  to  shoot  down  a  flying 
missile . 

Although  it  is  preferable,  it  is  not  necessary  to  cause 
the  ignition  of  the  warhead  in  order  to  destroy  a  flying 
target.  A  hole  of  the  dimensions  of  a  child's  fist  (3cm  x 
3cm  x  3cm)  is  generally  enough  to  cause  severe  damage  to  the 
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missile  consists  of  about  10  atoms,  so  that  about  3  MJ 
energy  is  needed  to  create  such  a  hole.  Limiting  the 
delivering  time  of  this  energy  to  2  to  3  seconds,  the  needed 
power  should  be  about  1  MW.  After  considering  diffraction, 
the  laser  beam  is  focused  at  the  target  in  an  area  of 

100  cm  so  that  the  intensity  is  about  10  kW/cm  .  In  order 
to  get  a  feeling  of  the  magnitude  of  this  intensity,  we  can 

mention  that  35  W/ cm  on  the  human  body  results  in  third 
degree  burns  with  destruction  of  the  dermis  and  epidermis. 
It  must  be  noted  that  delivering  time  cannot  be  diminished 
without  a  limit.  Smaller  delivery  times,  which  imply  very 
high  power,  may  cause  detrimental,  nonlinear  effects  such  as 
thermal  blooming.  These  effects  form  a  lens  in  the 
atmosphere  and  spread  the  laser  beam  diminishing  the  power 
at  the  target  dramatically.  Thermal  blooming  will  be 
discussed  in  the  following  section. 
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F.  PROPAGATION  OF  THE  LASER  BEAM  THROUGH  THE  ATMOSPHERE 

1 .  Transmittance 

The  atmospheric  components  (elements  and  compounds) 
strongly  affect  the  propagation  of  the  Laser  beam  and  can 
dramatically  diminish  the  required  optical  intensity  at  the 
target.  There  are  three  primary  processes  that  affect  the 
transmission  of  radiation  through  the  atmosphere: 
absorption,  scattering,  and  random  refractive  index 
fluctuations  (turbulence) . 

Absorption  is  a  quantum  process  where  an  atmospheric 
molecule  absorbs  the  energy  from  some  incident  photon  and 
occurs  when  certain  types  of  particles  interact  with  the 
propagating  electromagnetic  radiation.  Water,  vapor,  carbon 
dioxide,  ozone,  nitrous  oxide,  and  carbon  monoxide  are  the 
major  radiation  absorbers.  Scattering  is  the  redirection  of 
radiation  by  particles  in  the  air.  Photons  collide  with 
these  particles,  whereas  the  photon  energy  is  reradiated  in 
all  directions.  Turbulence  describes  the  time-space  varying 
temperature  and  pressure  inhomogenuites,  which  result  in 
considerable  index  of  refraction  fluctuations.  The  index  of 
refraction  variations  cause  the  direction  of  light 
propagation  to  bend  in  various  directions. 
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Absorption  and  scattering  are  grouped  together  under 
the  topic  of  extinction.  Extinction  is  defined  as  the 
attenuation  in  the  amount  of  radiation  passing  through  the 
atmosphere.  Absorption  and  scattering  are  described 
respectively  by  the  absorption  coefficient  s  and  scattering 
coefficient  p,  while  extinction  is  described  by  the 
extinction  coefficient  y  =  s  +  p  [  km  ] .  Transmittance  (Tr) 
through  the  atmosphere  over  a  distance  x  is  defined  as 

-\i-x 

Tr  =  e  ,  and  the  irradiance  at  a  target  at  distance  x  is 

I  (x)  =  IQe  p  =  I0Tr .  Both  £  and  p,  and  therefore  the 

transmittance  Tr,  strongly  depend  on  the  radiation 
wavelength. 

The  atmospheric  transmittance  over  a  wavelength  range 
is  extremely  crucial  to  directed  energy  because  it 
determines  the  choice  of  the  suitable  propagating 
wavelength.  Transmittance  depends  on  the  composition, 
density  and  pressure  of  the  atmospheric  components  together 
with  the  optical  line  strengths  and  line  widths  of  all  the 
spectroscopic  transitions  of  the  molecules  and  the  number, 
size  and  composition  distributions  of  the  particles. 

This  complex  problem  has  been  described  by  computer 
modeling  codes  such  as  LOWTRAN,  FASCODE  and  MODRAN.  MODRAN 
and  FASCODE  are  currently  the  most  widely  used  and  have  the 


35 


capability  to  handle  arbitrary  geometries  and  view 
conditions,  while  they  allow  the  user  to  select  from  several 
model  atmospheres,  aerosol  profiles,  and  other  weather 
conditions.  A  graph  of  a  typical  atmospheric  transmission 
in  a  maritime  aerosol  model  is  shown  in  Figure  15.  Note  the 
strong  dependence  of  the  transmittance  on  optical 
wavelength.  There  are  specific  windows  of  high 
transmittance  where  it  would  be  preferable  to  operate  a 
laser.  To  achieve  a  wavelength  in  one  of  these  windows,  a 
laser  must  be  tunable.  This  feature  is  unique  to  an  FEL; 
other  lasers  are  confined  to  specific  wavelengths. 


Figure  15.  Atmospheric  Transmittance.  (From  Ref.  5). 
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2. 


Nonlinear  Effects 


A  high  power  optical  beam  traveling  through  the 
atmosphere  is  also  subject  to  a  "nonlinear"  effect  where  the 
propagation  medium  is  modified  by  the  radiation.  This 
nonlinear  effect  is  called  thermal  blooming,  and  takes  place 
when  the  absorption  of  the  beam  energy  by  the  atmosphere 
causes  a  local  temperature  rise  followed  by  a  density 
change.  Since  the  refractive  index  depends  on  the 
atmosphere  density,  a  local  refractive  index  gradient  across 
the  laser  beam  cross-section  forms  a  lens  along  the  beam 
path.  This  effect  can  result  in  spreading,  distortion,  and 
bending  of  the  laser  beam,  while  it  generally  reduces  the 
power  density  at  the  target.  In  fact,  it  can  limit  the 
irradiance  to  a  value  independent  of  the  transmitted  power. 
Nonlinear  effects  depend  on  wavelength,  pulse  length, 
intensity,  and  such  medium  properties  as  composition  and 
density. 

As  a  result,  it  is  not  possible  to  increase  the  laser 
power  without  limit,  and  continue  to  have  the  laser  beam 
propagate  through  the  atmosphere.  Nonlinear  effects  can 
lead  to  dispersion  of  the  beam  power,  and  even  complete 
blockage . 
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6.  THE  FREE  ELECTRON  LASER  (FEL) 

Among  the  various  types  of  existing  lasers,  the  FEL 
turns  out  to  be  an  excellent  candidate  for  a  ship  self- 
defense  weapon  system.  It  can  be  easily  integrated  in  the 
ship's  structure  and  become  a  major  weapon  system.  The 
electrical  power  required  for  the  FEL's  operation  can  be 
provided  by  the  ship's  power  system  or  by  large  capacitors. 
The  weight  and  volume  of  all  the  FEL  components  together  are 

approximately  150  tons  and  140  m3  respectively  [Ref.  8], 
and  can  be  accommodated  in  a  ship  the  size  of  a  frigate  or  a 
destroyer. 

The  FEL  offers  significant  advantages  over  other  types 
of  lasers.  The  continuous  tunability  over  a  wide  range  of 
wavelengths  allows  the  opportunity  to  select  the  frequency 
that  best  propagates  through  the  atmosphere.  FELs  do  not 
produce  hazardous  and  toxic  wastes,  which  can  be  life 
threatening  for  the  ship's  crew  as  do  chemical  lasers. 
Finally  the  FEL  shows  the  possibility  of  high  wall  plug 
efficiency  ( ~ 1 0 % ) ,  making  good  use  of  the  ship's  valuable 
power  resources.  The  fact  that  it  operates  with  electrical 
power,  equips  the  FEL  with  an  infinite  magazine,  enabling  it 
to  minimize  reaction  time  and  engage  multiple  targets. 
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IV.  SIMULATION  RESULTS 


A.  THE  TAPERED  UNDULATOR 

At  normal  saturation,  the  electrons  lose  enough  energy 
to  drop  out  of  the  gain  spectrum  bandwidth.  At  resonance  we 

have  v  =  0  and  A  =  A0(l  +  K f  /  2y2  and  is  only  near  resonance 
that  the  coupling  between  the  electron  beam  and  the  light 
wave  is  significant.  If  A0  or  K  are  decreased  as  y 
decreases,  as  electrons  lose  energy,  resonance  would  be 
restored,  and  the  interaction  could  continue  to  higher 
optical  powers.  Increasing  A0  or  K  also  extends  the  FEL 
saturation  limit.  Decreasing  or  increasing,  tapering  or 
inverse  tapering,  the  undulator  wavelength  A0  or  K  along  z, 

modifies  the  simple  pendulum  equation  to  include  an 
additional  torque,  or  phase  acceleration  5,  so  that 

o  oo 

v  =  £  =  5  +  |a|  cos(£  +  cp) .  If  6>0,  we  have  the  case  of  the 


tapered" 

undulator, 

while 

for  5<0, 

the  undulator 

is 

inverse 

tapered" . 

The 

value  of 

5  is  given 

by 

5  =  -2nNAA0  /  A0  when  the  undulator  wavelength  is  decreased, 
or  by  5  =  -4nN  K2AB  /  B  (l  +  K2)  when  the  undulator  field 
strength  is  decreased. 
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Tapered  Undulator,  5=6n.  Figure  19.  Tapered  Undulator, 


Figures  16-19  show  the  computer-simulated  gain 
spectrum,  as  in  Figure  6,  of  an  FEL  with  dimensionless 
current  j=6  and  N=41  periods,  for  various  tapering  rates  5. 
We  can  see  that  as  the  tapering  rate  5  increases,  the  gain 
spectrum  is  no  longer  antisymmetric  about  resonance,  but 
distorts  as  the  bandwidth  broadens.  The  gain  in  weak 

optical  fields  ( a0  «  n )  is  bigger  in  the  untapered  case, 
while  in  the  strong  field  regime  the  gain  is  increased 
because  of  taper.  Figure  20  presents  the  peak  gain  of  an 
FEL  with  j=6  as  a  function  of  taper,  for  three  different 
optical  field  strengths  a0 .  In  weak  optical  fields,  a0  «  0 

when  the  magnitude  of  5  increases,  the  peak  gain  decreases. 
The  peak  gain  curve  is  symmetric  in  5  with  peak  value  at 
5=0.  When  the  optical  field  becomes  stronger,  the  shape  of 
the  curve  is  distorted;  the  peak  value  decreases  and  shifts 
to  negative  values  of  5.  Figure  20  shows  how  the  tapered 
undulators  have  higher  gain  in  strong  optical  fields  than 
the  untapered  case,  5=0.  In  strong  optical  fields  a0  =  40 
with  j  =  6,  the  maximum  gain  is  9%  achieved  with  5=-4n.  For 
5=0,  the  maximum  gain  in  strong  fields  is  only  4%.  For 
positive  values  of  5  and  a0  =40,  the  maximum  gain  of  8%  is 
achieved  with  5=10n. 
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Figure  20.  Plot  of  peak  gain  versus  tapering  rate. 

Figure  21  presents  a  plot  of  the  value  of  phase 
velocity  v0  giving  peak  gain  in  weak  fields  as  a  function 

of  the  tapering  rate  5.  In  weak  fields  and  for  moderate 
taper  rate  (0<5<10n) ,  the  gain  peaks  at  phase  velocity 
v0  =2. 4-5/2  .  Larger  taper  rates  and  strong  fields  alter 
this  relation. 
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Figure  21.  Position  of  Peak  Gain  versus  taper. 


B.  MULTIMODE  EVOLUTION  OF  SHORT  PULSES 

Multimode  behavior  is  simulated  by  introducing  a  short 
parabolic  electron  pulse  shape  into  the  optical  resonator 
and  following  the  evolution  of  the  optical  pulse  and  optical 
mode  as  a  function  of  the  number  of  round  trips  n,  through 
the  resonator.  Picosecond  long  electron  pulses  are  often 
used  in  the  FEL  oscillator  driven  by  an  RF  accelerator.  At 
resonance,  exactly  one  wavelength  of  light  A  passes  over 

an  electron  as  the  electron  passes  through  one  undulator 
wavelength  A0.  Over  the  whole  undulator  length,  N 
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wavelengths  of  light  pass  over  the  slower  electrons.  We 
define  the  "slippage  distance"  as  NA  and  is  the  distance 

that  the  electron  pulse  slips  behind  the  faster  optical 
pulse  as  they  travel  through  the  undulator.  It  is  the 
characteristic  length  over  which  electrons  and  light  can 
exchange  information  during  the  interaction  in  the 
undulator.  When  the  electron  pulse  length  is  comparable  to 
the  slippage  distance,  then  short  pulse  effects  are 
important  in  the  FEL  interaction. 

As  short  electron  pulses  enter  the  undulator,  short 
optical  pulses  start  from  spontaneous  emission  and  bounce 
between  the  resonator  mirrors  separated  by  a  distance  S, 
which  is  greater  than  the  undulator  length  L.  At  t=0  the 
electron  pulses  from  the  RF  accelerator  enter  the  undulator, 
while  the  rebounding  optical  pulses  arrive  to  overlap  the 
electron  in  time  intervals  of  2 S/c,  the  round-trip  bounce 
time  of  light  in  the  resonator.  We  define  "desynchronism" 
d  =  -2  AS/  (NA)  as  the  displacement  between  the  electron  and 
optical  pulses  at  the  beginning  of  the  undulator  (i=0) 
normalized  to  the  slippage  distance.  In  practice,  d  is 
adjusted  by  moving  one  of  the  resonator  mirrors,  thus 
altering  the  distance  S. 
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j=10  a2=l .  0  N=41  5=4n  Q=10  d=0 . 02 


-5  z  5  -47  V  47  0  n  2000 


Figure  22.  Short  Pulse  Evolution  with  5=4n  and  d=0.02. 

Figure  22  shows  a  sample  simulation  result  after  n=2000 
passes  with  taper  5=4n.  Several  parameters  are  printed 
across  the  top;  the  peak  current  j,  the  resonator  loss 
factor  Q,  the  pulse  width  az  ,  also  5,  and  N.  The  upper 
graphs  give  the  optical  field  shape  |a(z,  n)|  ,  the  optical 
power  spectrum  P(v,n),  and  the  electron  spectrum  f(v,n)  at 
the  final  pass.  The  shading  in  the  middle  graphs  show  how 
these  quantities  have  evolved  with  n.  The  scale  ranges  from 
0  in  black  to  the  maximum  value  shown  as  gray.  On  the 
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bottom  left,  the  longitudinal  profile  of  the  current 
density,  or  electron  pulse  j(z-i)  is  shown  for  reference  at 
dimensionless  times  i=0  (black)  and  x=l  (gray) .  The 
calculational  window  is  10  slippage  distances  long,  and 
electron  pulse  slips  behind  one  slippage  distance  in  that 
window.  The  z  parameter  in  this  Figure  is  dimensionless 
(  z/NA  )  .  The  bottom  center  graph  shows  the  weak-field  gain 
spectrum  for  reference,  and  the  right  bottom  graph  shows  the 
evolution  of  the  total  power  P(n).  The  power  is  the  square 
of  the  dimensionless  optical  pulse  amplitude  a  integrated 
over  the  pulse  length.  In  this  particular  case,  after  going 
through  some  transients,  the  FEL  reaches  a  steady  state 
within  approximately  n=100  passes.  Note  the  sudden  shift  of 
the  optical  power  spectrum  at  n  «  100;  as  n  increases,  the 
optical  field  strength  grows  and  distorts  the  gain  spectrum. 
As  seen  in  Figures  16-19,  the  gain  peaks  at  larger  values  of 
v  as  the  optical  field  a  becomes  large,  causing  the  power 
spectrum  to  shift  to  larger  values  of  v. 

Gain  in  weak  fields  is  described  by  G(t)  =  jv0i4  / 12  so 
that  there  is  no  gain  at  the  beginning  of  the  interaction. 
This  delay  in  gain  is  called  "lethargy".  At  later  times  the 
electron  pulse  becomes  bunched  and  the  light  pulse  is 
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distorted  because  gain  is  preferentially  deposited  on  the 
trailing  edge  of  the  pulse.  The  centroid  of  the  light  pulse 
travels  slower  than  c.  At  exact  synchronism  (d=0)  the 
electron  and  light  pulses  enter  the  undulator 
simultaneously.  Surprisingly,  the  steady-state  power  of  the 
FEL  in  this  case  is  0  because  the  light  pulse  drifts  away 
from  the  electron  pulse  over  many  passes.  Figure  23  shows 
the  pulse  evolution  of  such  a  case  where  the  optical  power 
evolves  to  zero  after  being  started  with  an  artificial  pulse 
of  light.  Introducing  a  small  value  of  desynchronism,  we 
compensate  for  the  slower  speed  of  the  light  pulse  and  we 
find  a  better  coupling  of  electron  and  optical  pulses. 
Considerable  power  is  achieved  in  the  steady  state,  as  shown 
in  Figure  22. 
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j=10  oz=l .  0  N=41  5=4n  Q=10  d=0 


-5  2  5  -47  V  47  0  n  2000 


Figure  23.  Short  Pulse  Evolution  with  5=4n  and  d=0. 

C .  TRAPPED-PARTICLE  INSTABILITY 

In  strong  optical  fields,  electrons  can  become  trapped 
in  deep  potential  wells  in  phase  space.  The  oscillation  of 
the  beam  current  that  is  trapped  in  these  wells  can  drive 
the  carrier  wave  unstable  and  cause  sideband  frequencies  to 
grow  from  noise.  As  the  FEL  reaches  high-power  saturation, 
the  height  of  the  separatrix  is  large  and  many  electrons  are 
trapped  in  the  closed  orbits  of  phase  space.  Trapped 
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electrons  execute  synchrotron  oscillation  with  the 
synchrotron  or  trapped-particle  oscillation  frequency 

vs  -  52j^ .  Sidebands  can  appear  around  the  fundamental 
at  v0  ±  vs ,  and  are  shifted  away  from  the  fundamental 
wavelength  by  AA  /  A  =  vs  /  2nN .  The  general  features  of 
trapped  particle  instability  depend  only  on  j  and  the  loss 
factor  Q.  Increasing  either  j  or  Q  can  result  in  the 
trapped  particle  instability.  The  trapped  particle 
instability  usually  occurs  for  small  desynchronism,  and  has 
mixed  effects  in  the  FEL.  Often,  the  experiment  would  like 
a  high  power  laser  that  has  a  narrow  spectrum  and  no 
sidebands.  But  the  presence  of  sidebands  means  that  high 
power  has  been  obtained  and  the  FEL  is  otherwise  operating 
well.  In  fact,  as  more  sidebands  develop,  even  more  power 
and  higher  efficiency  are  attained.  Figure  24  presents  the 
features  of  the  trapped  particle  instability.  The  optical 
pulse  (upper  left)  is  short  and  partially  modulated  with 
sharp  spikes.  A  small  sideband  is  also  evident  at  the 
optical  power  spectrum  P(v,n). 
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Figure  24.  Trapped  Particle  Instability. 


D.  SIMULATIONS  OF  THE  TJNAF  FEL  WITH  A  TAPERED  UNDUIATOR 

Experiments  using  the  TJNAF  FEL  have  explored  operation 
with  both  the  tapered  and  inversely  tapered  undulators  (Ref. 
15  and  16)  .  I  present  here  the  results  of  numerical 
simulations  using  the  TJNAF  experimental  parameters.  The 
TJNAF  FEL  has  N=41  undulator  periods  of  length  A0  =  2.7  cm 
each.  The  undulator  parameter  is  K=  0.98,  the  resonator 
loss  factor  is  Q=i0  (10%  loss  per  pass),  and  the  linear 


50 


taper 


rates 


are 


AK/K  =  5%  (5  =  4n)  ,  7.5%  (5  =  6n) , 


and  10%  (5  =  8n )  .  Multimode  simulations  describe  the 
evolution  of  short  pulses  in  the  far  infrared,  and  show  how 
positive  tapering  affects  single  pass  gain,  steady  state 
power,  and  the  electron  energy  spread  as  a  function  of 
desynchronism. 


1.  34.5  MeV  Electron  Beam 

In  this  case,  the  FEL  was  operated  with  short  0.5  ps 
electron  pulses  (length  le  =  150pm)  of  total  energy  E=34.5 

MeV,  with  a  0.25%  energy  spread.  The  peak  current  was 
1=50  A,  and  the  resulting  optical  wavelength  is  A  =  6pm. 
From  these  values,  the  dimensionless  parameters  where 
determined:  j=10  (dimensionless  current} ,  and 

c2  =  le  /  NA  =  1.0  (electron  pulse  length).  The  goal  was  to 

determine  the  steady  state  power,  the  steady  state  gain  in 
weak  fields,  as  well  as  the  electron  energy  spread  as  a 
function  of  both  5  and  d. 
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a)  Steady  State  Power 


Figure  25  summarizes  the  results  of  many 
simulations  obtained  using  numerous  values  of  desynchronism 
from  d=0  to  0=0.4,  for  taper  rates  5=  0,  +4n,  +6n,  and  +8n. 
The  number  of  passes  n  in  each  simulation  was  sufficient  for 
the  optical  power  to  evolve  to  steady  state.  The  steady- 
state  dimensionless  power  is  plotted  versus  desynchronism, 
for  different  tapering  rates  5,  giving  rise  to  the 


Figure  25.  Desynchronism  Curve  for  various  tapering  rates. 
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For  values  of  d  <  0.004,  the  laser  does  not  operate. 
But  at  slightly  large  value  of  d,  the  power  rises  sharply. 
For  small  d,  the  gain  is  small  so  that  a  large  number  of 
passes  n  is  required  to  achieve  steady-state  power.  For 
5=0,  the  sharp  peak  in  power  at  small  desynchronism  is 
accompanied  by  the  trapped-particie  instability.  This 
effect  is  reduced  as  taper  increases,  since  tapering  tends 
to  reduce  the  side  bands.  As  the  tapering  rate  is 
increased,  the  steady-state  power  decreases  and  the  peak 
power  of  each  curve  moves  slightly  to  larger  values  of 
de synchronism.  The  conventional  undulator  is  more  efficient 
than  the  tapered  undulators.  Furthermore,  higher  tapering 
rates  reduce  the  range  of  desynchronism  values  where  the 
laser  works.  For  no  tapering,  the  FEL  works  for  values  of 
desynchronism  as  big  as  d=0.38.  However,  a  tapering  rate  of 
5=8n  reduces  the  range  to  only  Ad=0.19.  For  5=0,  the  power 
peaks  at  d=0.004  with  the  trapped-particie  instability. 
This  is  expected  since  the  sidebands  contribute 
significantly  to  the  final  power.  For  5=8n,  the  steady- 
state  power  is  50%  smaller  and  peaks  at  d=0.015.  For  larger 
d,  the  power  diminishes  significantly  as  the  taper 
increases.  The  flat  sections  of  the  desynchronism  curves 
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since  they  are 


are  desirable  regions  to  operate  the  FSL, 
predictable  and  stable.  Limit-cycle  behavior  is  observable 
for  5=0  between  d=0.02  and  d=0.07,  and  for  6=8n  between 
d=0.025  and  d=0.03.  For  these  regions,  shown  on  the  graph 
as  large  circles,  only  the  peak  power  of  the  steady  state 
power  is  shown.  Limit-cycle  behavior  causes  oscillations 
of  the  power  P(n),  the  power  spectrum.  P(v,n),  and  the 
electron  spectrum  f(v,n).  The  steady-state  power  oscillates 
periodically  by  as  high  as  50%,  between  the  peak  and  the  low 
values.  Limit-cycle  behavior  occurs  when  trapped  particles 
in  strong  fields  combine  with  short  optical  pulses.  The 
modulation  caused  by  the  oscillation  of  the  trapped  current 
continually  modifies  the  shape  of  the  short  optical  pulse. 
The  different  pulse  shapes  have  different  powers  ?  and 
spectra  P(v),  causing  oscillations  as  subpulse  structures 
march  through  the  optical  pulse  envelope.  Figure  26  shows 
an  example  of  limit  cycle  behavior  for  5=8rr  and  d=0.026. 
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j=10  az=l .  0  N=41  6=4n  Q=10  d=0 . 12 


-2  z  2  -50  V  50  0  a  2000 


Figure  26.  Limit-cycle  behavior. 


b)  Steady  State  Gain 

In  order  to  determine  the  steady-state,  weak  field 
gain  as  a  function  of  desynchronism,  the  same  simulations 
were  used.  Instead  of  plotting  power  as  a  function  of  n  in 
the  lower  right  graph  of  each  simulation  output,  P(n)  is 
replaced  by  gain  G(n),  as  shown  in  Figure  27.  The 
simulation  was  run  only  long  enough  for  gain  to  reach  a 
steady-state  value. 
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j=10  o2= 1.0  N=41  6=4n  Q=oo  d=0 . 12 


-5  2  5  -47  V  47  -47  V  47 


-5  2  5  -47  V  47  0  n  30 


Figure  27.  FEL  Short  Pulse  arid  Gain  Evolution. 

Since  we  are  interested  in  the  weak  field  gain, 
the  optical  field  a  must  be  kept  less  than  unity  by  starting 

with  an  initial  optical  field  ac  *  ICf2 .  The  loss  factor  is 

Q  =  10i0  »  oo ,  so  that  the  mirror  losses  (1/Q)  would  be 
negligible.  The  power  and  electron  spectra  are  narrow, 
since  the  optical  field  is  small. 
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Figure  28.  Steady  State  Gain  versus  desynchronism 


In  Figure  28,  the  weak-field,  steady-state  gain 
is  plotted  as  a  function  of  desynchronism  d,  for  the  given 


tapering 

rates 

5. 

Clearly  the  gain  curves 

are 

much 

different 

from 

the 

desynchronism  curves  and 

peak 

for 

moderate 

values 

of  d. 

When  tapering  increases. 

the 

peak 

gain  of  each  curve  shifts  to  smaller  values  of  d.  For  no 
tapering  (5=0),  the  peak  gain  is  85%,  at  d=0.14.  When  5=8n, 
the  peak  gain  drops  to  only  12%,  at  d=0.03.  The  most 
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common  feature  of  the  desynchronism  curve  is  that  the 
operating  range,  (Ad),  decreases  as  5  increases.  Note  that 
when  the  gain  drops  to  Q=10%,  the  final  power  drops  to  0 
because  the  FEL  loss  per  pass  exceeds  the  gain  per  pass. 

c)  Electron  Energy  Spread 

When  the  power  reaches  steady  state,  we  measure 
the  full  width  Av  of  the  final  electron  spectrum,  which  is 
the  upper  right  graph  of  the  computer  simulation  output  (see 
Figure  24) .  From  Av  we  can  determine  the  fractional 
electron  energy  spread  using  Ay/y  =  Av  /  4nN  where  N=41.  In 
the  TJNAF  EEL  design,  the  electron  beam  exiting  the 
undulator  is  redirected  by  bending  magnets  and  fed  back  into 
the  accelerator.  This  recirculation  allows  energy  recovery 
from  the  electron  beam,  and  significantly  increases  the  FEL 
efficiency.  The  induced  electron  energy  spread  is  crucial 
to  this  process,  because  beam  recirculation  is  feasible  only 
if  the  full  electron  energy  spread  is  less  than  6%.  In 
Figure  29,  the  fractional  electron  spread  induced  by  the  FEL 
interaction  is  plotted  as  a  function  of  d  and  5.  The  curves 
appear  much  like  the  power  curves  in  Figure  25  with  initial 
sharp  peaks  of  7%-8%  near  d=0.01.  However  for  d  >0.025,  all 
the  energy  spreads  are  less  than  6%. 
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Figure  29.  Electron  Energy  Spread. 


Between  d=0  and  d=0.05,  where  we  have  the  maximum 
output  power,  the  induced  energy  spread  is  significantly 
smaller  for  the  conventional  undulator  than  the  tapered  one. 
For  d=0.025,  the  energy  spread  was  4.5%  for  the  untapered 
undulator  and  around  6%  for  all  the  other  tapers.  The 
tapered  undulator  does  not  appear  to  be  desirable  for 
recirculating  the  electron  beam.  After  d=0.05,  the  energy 
spread  drops  to  less  than  4%  for  all  the  tapering  rates. 
With  short  optical  pulse,  the  head  and  tail  of  the  pulse 
have  weak  optical  fields,  so  that  tapering  is  not  optimum 
along  the  pulse.  This  is  the  reason  the  energy  spread  is 
greater  in  the  tapered  undulators. 
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2. 


47.5  MeV  Electron  Beam 


defined: 


a) 

Steady- State  Power 

In 

this 

case,  the 

dimensionless 

parameters  are 

j=7 

and 

oz  =  1.8  . 

The  optical 

wavelength  was 

A  =  3iam.  The  desynchronism  curves  were  again  determined  as 


shown  in  Figure  30.  The  curves  are  very  similar  with  those 
using  the  32.5  MeV  (Fig.  25).  However,  the  operating  ranges 
and  the  final  power  are  smaller  because  of  the  smaller  value 
of  j.  For  taper  5=8n,  the  laser  barely  works. 


Figure  30.  Desynchronism  Curve  for  47.5  Mev  electron 
beam  and  various  tapering  rates. 
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The  operating  width  of  the  desynchronism  curves, 
Ad,  as  a  function  of  taper  has  been  measured  in  experiments 
conducted  with  the  IR  Demo  FEL  at  Jefferson  Lab  [Ref  16]  . 
The  FEL  parameters  used  in  the  simulations  were  the  same  as 
those  of  the  experiments  (A  =  3pm,  N=41,  Q=10,  j=7  and 

oz  =  1.8).  In  Figure  31,  the  relative  width  of  the 
desynchronism  curve  is  plotted  as  a  function  of  taper,  using 
both  experimental  and  computer  simulation  data.  The 
relative  width  is  the  maximum  value  of  desynchronism  where 
the  FEL  operates  for  each  taper,  divided  by  the  width  of  the 
desynchronism  curve  of  the  untapered  undulator  (5=0) .  Both 
simulations  and  experiments  show  that  the  operating  range 
decreases  as  the  taper  increases.  There  is  good  agreement 
between  theory  and  experiment. 


Figure  31.  Relative  width  of  desyncronism  curve  versus  taper  5. 
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The  trapped-particle  instability  and  limit-cycle 
behavior  were  again  characteristic  features  of  the  simulation 
results.  Figure  32  presents  the  simulation  output  with  j=7 
and  az  =  1.8  (47.5  MeV  Energy),  5=0,  and  d=0.025.  The 


simulation  was  run  only  for  n=300  passes  in  order  to  better 
present  the  details.  The  oscillation  of  the  total  power, 
the  optical  pulse  distortion  and  the  sideband  in  the  optical 


power  spectrum,  are  clearly  shown. 
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Limit-cycle  behavior  in  the  final  power  was 
observed  for  5=0,  4n,  8n  and  not  for  5  =  6n.  Figures 
33-35  show  the  desyncronism  curves  of  Figure  30  in  order  to 
more  clearly  presentthe  details  of  limit-cycle  behavior. 


Figure  33.  Desynchronism  Curve  for  6=0. 
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Figure  34.  Desynchronism  Curve  for  5=4n. 


The  upper  right  graphs  of  Figures  33  and  34  show 
blown  up  section  of  the  desynchronism  curve  section  where 
the  limit-cycle  behavior  occurs.  Since  there  are  periodic 
oscillations  of  the  total  power,  the  maximum  and  minimum 
power  is  indicated  in  the  graph.  The  mean  power  is  the 
average  of  the  maximum  and  the  minimum.  The  desynchronism 
curves  in  Figure  30  plotted  the  mean  power.  Figure  35  shows 
the  desynchronism  curve  for  5=8n,  emphasizing  the  limit- 
cycle  behavior. 
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Figure  35.  Desynchronism  Curve  for  5=8n. 

b)  Steady-State  Gain 

The  same  procedure  as  the  one  followed  in 
paragraph  lb  was  followed  in  order  to  determine  the  steady- 
state,  weak-field  gain  as  a  function  of  taper  and 
desynchronism.  The  results  are  presented  in  Figure  36  and 
are  similar  to  those  with  the  32.5  Mev  beam  energy 
(Figure  29)  .  Still  the  operating  ranges,  as  well  as  the 
gains,  are  smaller. 
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Figure  36.  Steady  State  Gain  versus  desynchronism  for  47.5 
Mev  electron  beam. 


c)  Electron  Energy  Spread 

Figure  37  shows  the  fractional  electron  spread  as 
a  function  of  d  and  5.  The  results  are  again  similar  to 
those  with  32.5  MeV  electron  beam  energy.  The  curves  have 
the  same  trend  as  the  power  curves  (Fig.  30)  with  a  peak 
near  d=0.01.  The  peak  fractional  energy  spread  is 
approximately  1%  smaller  than  for  the  34.5  MeV  beam  because 
of  the  smaller  value  of  j . 


66 


Figure  37.  Electron  Energy  Spread  for  42.5  MeV  Electron  Beam. 

Initial  peaks  of  7%  induced  energy  spread  are 
evident  at  the  start  up  region  (d=0.007).  For  d  >  0.015, 
the  energy  spread  drops  to  less  than  6%,  implying  the 
possibility  of  beam  recirculation.  In  the  region  of  high 
steady-state  power (0.01<  d  <0.04)  the  conventional  undulator 
(5=0)  has  almost  1%  less  energy  spread  than  the  tapered 
(5=4n,6n).  Thus,  the  untapered  undulator  is  a  better 
candidate  for  beam  recirculation.  The  short  optical  pulses, 
with  weak  fields  at  the  head  and  tail,  cannot  be  easily 
designed  to  reduce  energy  spread. 


E.  CONCLUSIONS 

Experiments  were  conducted  with  the  TJNAF  IR  demo  FEL 
(Ref.  16)  in  order  to  test  its  operation  with  various  taper 
rates  5,  and  desynchronism  values  using  the  following 


parameters : 


Undulator 

Periods 

N 

41 

41 

Wavelength 

A 

3  pm 

6  pm 

Electron  Beam 
Energy 

E 

47.5  MeV 

34.5  MeV 

Table  3.  Parameters  used  in  the  FEL  experiments. 


In  this  thesis,  computer  simulations  were  used  in  order 
to  explore  the  operation  of  the  TJNAF  FEL  with  the 
conventional  and  several  tapered  undulators  using  the  same 
parameters  as  in  experiments.  The  corresponding 
dimensionless  parameters  used  in  the  simulations  are  listed 


in  Table  4. 


Electron  Beam 
Energy 

E 

47.5  MeV 

34.5  MeV 

Dimensionless 

Current 

j 

7 

10 

Resonator 
loss  factor 

Q 

10 

10 

Pulse  length 

oz 

1.8 

1.0 

Table  4.  Parameters  used  in  the  simulations. 


The  desynchronism  curves  were  determined  at  each 
electron  beam  energy,  and  present  a  sharp  peak  at  small 
values  of  desynchronism  accompanied  by  the  trapped-particle 
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instability.  Trapped-particle  instability  was  reduced  for 
larger  values  of  taper,  since  tapering  suppresses  the 
sidebands.  Limit  cycle  behavior  was  also  observed  for  most 
tapering  rates.  Although  the  general  shape  of  the  curves  is 
very  similar  for  each  beam  energy,  the  final  power  and  gain 
were  significantly  smaller  for  the  47.5  MeV  electron  beam 
than  for  the  34.5  MeV  electron  beam,  because  of  the  lower 
value  of  j . 

The  tapered  undulator  turned  out  to  be  less  efficient 
than  the  conventional  undulator,  demonstrating  reduced  gain 
and  power,  and  increased  energy  spread  to  the  electron  beam. 

The  width  of  the  desynchronism  curve  decreased  as  the 
tapering  rate  is  increased,  in  agreement  with  the 
experiments.  A  publication  is  being  prepared  comparing  the 
experimental  (Ref.  16),  and  simulation  results. 
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V.  EXPERIMENTAL  RESULTS  OF  DAMAGE  INDUCED  TO  MATERIALS  BY 

THE  TJNAF  FEL 

A.  LASER  MATTER  INTERACTION 

The  interaction  between  lasers  and  matter  is  a 
complicated  issue  generally  described  by  non-linear  and 
unpredictable  effects.  The  laser  beam  has  the  unique 
ability  to  deliver  very  high  power  per  unit  area.  When  high 
power  laser  radiation  falls  on  a  target,  the  part  of  the 
beam  that  is  absorbed  begins  to  heat  the  target  surface  very 
rapidly  to  its  melting  temperature.  This  melting  process 
then  penetrates  progressively  into  the  material.  Many 
physical  processes  govern  the  damage  caused  to  the  material 
by  the  laser  power  including  power  absorption,  power 
reflection,  heat  conduction,  and  heat  diffusion. 
Furthermore,  a  large  number  of  parameters  play  a  major  role 
in  these  processes  such  as  material  density  and  heat 
capacity,  as  well  as  the  irradiation  wavelength,  power 
density,  peak  power,  and  possibly  pulse  characteristics. 

A  good  knowledge  of  these  mechanisms  helps  one 
understand  the  capabilities  and  limitations  of  the  directed 
energy,  and  allows  better  control  of  the  damage  caused  by 
the  laser.  Controlled  damage  has  many  industrial 
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applications  such  as  the  creation  of  thin  coatings, 
electronic  component  fabrication,  precise  drilling,  cutting, 
etc.  However,  when  using  the  high-power  laser  as  a  weapon 
to  shoot  down  missiles,  precision  and  symmetry  of  the 
damaged  area  are  not  the  issue.  Our  goal  is  to  cause  the 
maximum  possible  damage  on  a  rapidly  moving  missile;  this 
requires  precise  tracking  of  the  laser  beam  through  a 
turbulent  atmosphere 

B .  SCALING 

As  stated  in  chapter  III,  an  intensity  of  10  kW  /  cm2 
over  a  6  cm  radius  spot  for  a  few  seconds  dwell  time  is 
needed  to  burn  an  adequate  hole  in  a  missile.  In  these 
damage  experiments,  the  TJNAF  FEL  used  only  several  hundred 
watts.  In  order  to  achieve  the  desired  intensity  of  10 

kW  /  cm  ,  the  laser  beam  has  to  be  focused  with  a  lens  to  a 
much  smaller  spot  size.  For  average  power  P  and  spot  area 
A,  the  resultant  intensity  is  I  =  P  /  A  .  Given  a  specific 
power,  we  can  adjust  the  spot  size  to  achieve  the  desired 
intensity.  Studying  a  spot  radius  on  the  order  of  0.1  cm 

and  using  500  W  power,  we  can  achieve  10  kW  /  cm2  intensity. 
The  goal  of  this  scaling  is  to  develop  guidelines  that  will 
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reliably  predict  the  damage  caused  by  a  high  power  laser 
over  a  large  area. 

Nevertheless,  the  spot  size  cannot  be  arbitrarily 
small,  or  the  scaling  will  not  work.  Each  material  is 
characterized  by  a  parameter  called  the  thermal  diffusion 
length  D,  which  represents  the  distance  required  for  the 
temperature  to  drop  by  a  factor  of  1/e,  and  determines  the 
ability  of  the  material  to  absorb  and  transport  heat.  If  D 
is  greater  than  the  laser  beam  spot  size  diameter,  the  heat 
deposited  by  the  laser  beam  will  diffuse  away  in  less  time 
than  it  takes  the  material  to  melt.  In  order  to  obtain 
effective  heating,  and  cause  melting  of  the  material 
irradiated,  the  thermal  diffusion  length  should  be  smaller 
than  the  beam  diameter.  Also,  the  damage  hole  cannot  be 
much  greater  than  the  laser  beam  radius.  Material  cannot 
easily  escape  a  deeply  damaged  hole,  and  the  damage 
mechanism  may  be  complicated  by  the  material  ejecting  the 
hole  and  flowing  through  the  beam. 


C .  EXPERIMENTAL  PROCEDURE 

The  sample  materials  used  in  the  experiments  were 
Polyimide  Fiberglass  and  F2  Epoxy.  The  samples  provided  by 
the  Naval  Research  Laboratory  are  the  same  as  those  used  in 
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similar  experiments  conducted  in  March  1999.  The  results  of 
the  March  1999  irradiations,  which  involved  three  holes  in 
each  sample  labeled  1,2,  and  3,  were  presented  in  Reference 
[14],  and  some  of  them  will  be  mentioned  in  this  thesis. 
The  results  presented  by  this  thesis  were  obtained  by  two 
different  experiments  conducted  in  August  1999  and  March 
2000.  The  results  from  both  experiments  will  be  presented 
together  in  order  to  make  the  analysis  more  comprehensive. 

1.  August  1999  Experiment 

Each  sample  was  irradiated  12  times  (4  sets  of  3 
irradiations  with  the  same  parameters)  with  a  FEL  beam  of 
wavelength  A=3.1  pm,  pulse  repetition  frequency  18.7  MHz  and 
average  power  P=100W±5W.  The  Polyimide  Fiberglass  was 
irradiated  first,  followed  by  the  F2-Epoxy.  Measurements 
were  made  with  samples  placed  downstream  of  a  calcium 
fluoride  lens  with  a  measured  back  focal  length  of  137.6  mm 
for  3  pm  wavelength.  A  camera  was  set  up  to  observe  the 
front  and  back  surface  of  the  samples.  Two  sets  of  three 
irradiations  were  first  made.  The  average  intensity  was  500 

W  /  cm  ,  which  was  achieved  by  focusing  the  beam  to  a  spot 
of  0.25  cm  radius.  Three  identical  irradiations  were  made 
initially  with  no  airflow.  Then  three  more  irradiations  were 


74 


made  with  a  wind  speed  of  60  mph  across  the  front  face  of 
the  samples.  Then,  adjusting  the  beam  radius  to  0.087  cm, 

the  intensity  was  set  to  1=10  kW  /  cm  and  the  same  set  of 
measurements  were  repeated. 

2 .  March  2000  Experiment 

The  goal  of  these  measurements  was  to  maintain  the 

average  intensity  of  10  kW  /  cm  ,  but  using  higher  laser 
power  with  a  larger  spot  size.  Each  sample  was  again 
irradiated  6  times  (2  sets  of  3  irradiations)  with  a  FEL 
beam  of  wavelength  A=3.1  pm,  pulse  repetition  frequency 
37.425  MHz  and  average  power  P=500W±10W.  The  Polyimide 
Fiberglass  was  irradiated  first,  followed  by  the  F2-Epoxy. 
Measurements  were  made  with  samples  placed  downstream  of  a 
calcium  fluoride  lens  with  a  measured  focal  length  of  235.7 
mm  for  3  pm  wavelength.  A  camera  was  again  set  up  to 
observe  the  front  and  back  surface  of  the  samples. 

Irradiations  were  made  with  1=10  kW  /  cm  average  intensity, 
which  was  achieved  by  focusing  the  beam  to  a  spot  of  0.13  cm 
radius.  Three  identical  irradiations  were  made  with  no 
airflow,  and  then  three  more  irradiations  with  a  wind  speed 
of  85  mph  across  the  front  face  of  the  samples. 
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D .  RESULTS 

1 .  Sample  #1  -  Polyimide  Fiberglass 

The  sample  had  dimensions  11.4  cm  by  10.1  cm  with  2  mm 
thickness.  In  Figure  38,  we  see  a  photo  of  the  front  side 
of  the  Polyimide  Fiberglass  sample  after  all  sets  of 
irradiations.  Irradiations  labeled  by  numbers  1,2  and  3 
were  conducted  in  March  1999.  Irradiations  labeled  4  to  15 
were  conducted  in  August  1999,  and  16  to  21  were  conducted 
in  March  2000.  All  irradiations  were  done  three  times  with 
the  same  parameters  in  order  to  get  more  accurate 
measurements.  The  actual  results  came  from  the  mean  value 
of  the  three  measurements.  Irradiations  7,8,9,10,11,12,16 
17  and  18  were  done  with  no  airflow,  while  in 
4,5,6,13,14,15,19,20,  and  21  there  was  wind  present.  In 
Figure  39,  we  see  a  photo  of  the  backside  of  the  Polyimide 
Fiberglass  sample  showing  that  all  of  the  irradiations  with 

intensity  of  10  kW  / cm2  completely  penetrated  the  sample. 
On  the  other  hand,  the  irradiations  with  an  intensity  of 

500  W  /  cm2  (irradiations  10  to  15  of  Figure  38)  did  not 
penetrate  the  sample.  Irradiation  results  are  summarized  in 
Table  5.  Each  row  represents  a  set  of  three  irradiations. 
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The  hole  diameters  have  variation  of  15%  in  each  set  of 


three  holes.  This  is  due  to  the  slight  variation  of  the 
exposure  time,  which  is  on  the  order  of  1-2  seconds. 

The  presence  of  airflow  resulted  in  «15%  bigger  hole 
diameter  and  penetration  depth  rate,  than  those  attained 
without  airflow.  The  irradiations  conducted  with  500  W 
average  power  (runs  16  to  21),  caused  3.5  times  higher 
penetration  depth  rate  than  the  ones  conducted  with  100  W, 

and  the  same  intensity  ( kW  /  cm2 )  .  This  is  probably  because 
the  higher  power  allowed  larger  spot,  so  that  damaged 
material  could  more  easily  escape  the  larger  damage  hole. 

The  damage  pattern  of  runs  16  to  21  is  not  circular  as 
expected  but  somewhat  elliptical,  which  is  very  clear  in 
Figures  44,  45,46  and  47.  The  following  figures  present  a 
closer  caption  of  the  damage  caused  by  the  irradiations. 


79 


Figure  40.  Poiyimide  hole  12. 

2 

Intensity  500  W  /  err;  ,  no  air. 
Average  Power  100  W. 

Spot  Radius  2.5  mm. 


Figure  41.  Poiyimide  hole  15. 

2 

Intensity  500  W  /  cm  with  air. 
Average  Power  100  W. 

Spot  Radius  2.5  mm. 


Figure  42.  Poiyimide  hole  8. 

2 

Intensity  10  kW  /  cm  ,  no  air. 
Average  Power  100  W . 

Spot  Radius  0.87  mm. 


Figure  43.  Poiyimide  hole  5. 

2 

Intensity  10  kW  /  cm  ,  with  air. 
Average  Power  100  W. 

Soot  Radius  0.87  mm. 
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Figure  44.  Polyixrd.de  bole  18. 

2 

Intensity  10  kW  /  cm  ,  no  air. 
Average  Power  500  W. 

Spot  Radius  1.25  mm. 


Figure  45.  Polyimide  hole  20. 

2 

Intensity  10  kW  /  cm  ,  with  air. 
Average  Power  500  W. 

Spot  Radius  1.25  mm. 


2  mm 


Figure  46.  Polyimide  exit  -hole  18. 

2 

Intensity  10  kW  /  cm  ,  no  air. 
Average  Power  500  W. 

Spot  Radius  1.25  mm. 


Figure  47.  Polyimide  exit  hole  20. 

2 

Intensity  10  kW  /  cm  ,  with  air. 
Average  Power  500  W. 

Spot  Radius  1.25  mm. 
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In 

all  cases  a  raised 

lip  of  melted 

material 

is 

observed 

around  the  face  of 

the  entrance 

hole . 

The 

dimensions  of  the  lip  are  approximately  0.2  mm  height  and 
1  mm  width.  However,  when  airflow  is  present  the  lip  tends 
to  be  smaller,  possibly  because  the  airflow  removes  the 
debris,  and  the  evaporated  material  that  actually 
contributes  to  the  formation  of  the  lip.  That  would  explain 
the  bigger  diameter  of  the  holes  made  in  the  presence  of 
wind,  since  the  diameter  is  measured  from  the  inside  part  of 
the  lip. 

The  charred  region  extends  approximately  2  mm  around 
the  lip  when  there  is  no  wind.  With  wind  this  area  is  much 
smaller.  The  charred  region  around  the  exit  hole  extends  to 
1  to  1.5  mm  in  all  cases,  as  the  backside  of  the  sample  was 
not  exposed  to  the  airflow.  As  exposure  time  increased  the 
radial  extent  of  the  damaged  area  increased  and  more  melted 
material  was  deposited  around  the  hole.  After  the 
irradiation  stopped,  there  was  a  period  of  almost  3  sec  that 
the  material  is  still  hot  and  melting.  The  wind  tends  to 
cool  down  the  material,  decreasing  this  time  almost  in  half 
and  resulting  in  less  melted  material. 

Investigation  with  a  microscope  reveals  that  there  is 
no  evidence  of  melted  or  rehardened  material  inside  the 
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holes . 


It  is  also  evident  that  the  damage  is  most 


significant  in  the  center  of  the  hole,  diminishing  radially 
outward,  which  supports  expectation  of  a  Gaussian  shape  of 
the  laser  beam  intensity. 

The  following  table  presents  the  results  of  the 
experiment  conducted  in  March  1999  (Holes  1,2,3). 


Average 

Power 

Average 

Intensity 

PRF 

Spot 

radius 

Wind 

Penetration 

Depth 

(W) 

(kW/cm2) 

(MHz) 

(pm) 

(mm) 

(mph) 

Rate 

(mm/sec) 

100 

10 

0.87 

NO 

0.9 

Table  6.  Irradiation  results  of  March  1999  experiment  on 
Polyimide  Fiberglass. 

The  above  results  have  the  same  irradiation  parameters 
with  those  of  runs  7,8,9  of  Table  5  except  for  the  higher 
Pulse  Repetition  Frequency  (PRF)  of  37.4  MHz  and  the  longer 
wavelength  A=4 . 8  pm.  It  appears  that  the  shorter  wavelength 
in  Table  5  (A=3.1  pm)  combined  with  the  lower  PRF  is  much 
more  effective,  resulting  in  60%  higher  penetration  depth 
rate  than  in  table  6.  The  lower  PRF  apparently  causes  more 
damage,  since  it  results  in  more  energy  per  pulse  for  the 
same  average  power. 
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2 .  Sample  #2  -  F2  Epoxy 

The  sample  had  dimensions  11.5  cm  by  10  cm  with  1.5  mm 
thickness,  attached  to  a  1.6  cm  thick  polyethane  foam 
backing  which  is  clearly  seen  in  Figure  49.  In  Figure  48, 
we  see  a  photo  of  the  front  of  the  sample  after  all  sets  of 
irradiations.  Irradiations  labeled  1,  2  and  3  were  those 

conducted  in  March  1999.  Irradiations  4  to  15  were 
conducted  in  August  1999,  and  16  to  21  were  the  last  ones 
conducted  in  March  2000.  Following  the  same  procedure  as 
with  polyimide,  all  the  irradiations  were  done  three  times 
with  the  same  parameters  in  order  to  get  more  accurate 
measurements.  The  results  presented  here  came  from  the  mean 
value  of  the  three  measurements.  Irradiations  7,8,9,10, 
11,12,16,17,  and  18  (Fig.  48)  were  done  with  no  airflow 
while  in  4,5,6,13,14,15,19,20,21  there  was  wind  present.  In 
Figure  50,  we  see  a  photo  of  the  backside  of  the  sample 

showing  that  all  of  the  10  kW  /  cm2  irradiations  completely 

penetrated  the  sample.  On  the  other  hand  the  500  W  /  cm2 
intensity  (runs  4  to  9)  did  not  penetrate,  but  caused  more 
extensive  surface  damage  due  to  the  bigger  spot  radius. 
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Figure  48.  F2  Epoxy  (front  view). 


Figure  49.  F2  Epoxy  (side  view) . 
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Figure  50.  F2  Epoxy  (back  view). 


The  presence  of  the  foam  layer  at  the  backside  of  the 
sample  made  the  measurements  of  the  exit  holes  diameter 
unreliable.  A  closer  caption  of  the  damage  is  presented  in 
Figures  51-56. 

Irradiation  results  are  summarized  in  Table  7.  Each 
row  of  the  table  represents  a  set  of  three  irradiations, 
conducted  with  the  same  characteristics. 
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Figure  51.  Epoxy,  hole  7. 

2 

Intensity  500  w  /  c~  ,  no  air. 

Average  Power  IOC  W .  Spot  Radius  2.5  mra. 


Figure  52.  Epoxy,  hole  5. 

2 

Intensity  500  w  /  err.  ,  with  airflow. 
Average  Power  100  W.  Spot  Radius  2.5  mn. 
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Figure  53.  Epoxy,  hole  16. 

Intensity  10  kW/cm2 ,  no  air. 

Average  Power  500  W.  Spot  Radius  1.25  mm. 


Figure  54.  Epoxy,  hole  20. 

Intensity  10  kW/cm2,  with  airflow. 

Average  Power  500  W.  Spot  Radius  1.25  mm. 
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Figure  55.  Epoxy,  hole  12. 

Intensity  10  kw/cra2,  no  air. 

Average  Power  100  W.  Spot  Radius  0.87  nan. 


Figure  56.  Epoxy,  hole  15. 

Intensity  10  kW/crn2,  with  airflow. 
Average  Power  100  W .  Spot  Radius  0.87  xrnn. 
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The  hole  diameters  have  variation  of  15%  in  each  set  of 


three  holes.  This  is  again  due  to  the  slight  variance  of 
the  exposure  time  of  the  irradiations,  which  is  on  the  order 
of  1  to  2  seconds.  The  presence  of  airflow  did  not  seem  to 
increase  the  extent  of  the  damage.  The  entry  hole  diameter 
was  actually  decreased  by  10%  to  30%. 

The  average  power  of  500  W  (runs  16  to  21),  resulted  in 
2.5  times  higher  penetration  rate  than  the  100  W  power  (runs 
10  to  15) ,  and  caused  a  slightly  elliptical  damage  pattern. 
During  the  irradiation,  it  was  observed  that  flames,  smoke, 
and  debris  covered  the  entrance  hole.  After  the  irradiation 
stopped  the  material  was  still  burning  for  almost  3  seconds, 
which  caused  extra  charring  and  melting  of  the  sample.  When 
airflow  was  applied,  the  time  decreased  by  half.  The 
charred  region  extends  approximately  0 . 5  mm  around  the 
entrance  hole  with  wind  present  and  1  mm  without  wind. 
Examination  of  the  holes  with  a  microscope  revealed  more 
roughness  than  the  polyimide.  This  was  probably  caused  by 
deposited  debris  and  charred  material.  The  Gaussian  beam 
caused  the  same  damage  pattern  with  the  polyimide  sample, 
being  more  intensive  at  the  center  of  the  spot.  However  the 
penetration  rates  were  observed  to  be  2  to  3  times  smaller. 
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Table  8  presents  the  results  of  the  experiment 
conducted  in  March  1999  for  the  Epoxy  sample  (Holes  1,2,3). 


Average 

Power 

Average 

Intensity 

PRF 

A 

Spot 

radi 

Wind 

Penetration 

Depth 

(W) 

(kW/cm2 

(MHz) 

(pm) 

us 

(mm) 

(mph) 

Rate 

(mm/sec) 

100 

10 

37.4 

4.825 

0.87 

NO 

0.1 

Table  8.  Irradiation  results  of  March  1999  experiment  on  Epoxy 
Comparing  these  results  with  runs  10  to  12  of  Table  7, 
it  is  clear  again  that  the  A=3.1  pm  wavelength  combined  with 
the  lower  Pulse  Repetition  Frequency  was  more  effective,  and 
resulted  in  6  times  higher  penetration  depth  rate. 


E.  CONLUSIONS  -  SUGGESTIONS 

When  higher  laser  power  is  available  in  the  future,  the 
same  irradiations  should  be  conducted  with  even  larger  spot 
radii  in  order  to  compare  the  results  and  establish  scaling 
laws.  It  is  suggested  that  thicker  samples  should  also  be 
tested  in  order  to  determine  if  the  penetration  rate  stays 
the  same  after  the  laser  beam  has  penetrated  the  material  a 
few  mm. 
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Airflow  did  not  have  a  significant  effect,  possibly 
because  of  its  relatively  low  speed.  It  is  suggested  that 
in  future  experiments  the  wind  speed  should  be  much  higher 
in  order  to  simulate  more  realistically  the  conditions  of  a 
missile  flying  with  1  to  2  Mach  speed  (1000  mph) . 

Furthermore,  a  way  to  make  the  measurements  of  burn 
through  times  more  accurate  should  be  established,  as  they 
are  important  in  determining  penetration  depth  rates. 

From  the  analysis  of  the  results,  it  is  apparent  that 
the  penetration  depth  rate  is  proportional  to  the  laser 
intensity.  When  the  intensity  is  changed  by  a  factor  of  20 

(from  500  W  /  cm2  to  10  kW  /  cm2 )  the  penetration  depth  rate 
is  also  changed  by  almost  the  same  factor  in  both  samples. 

The  A=3.1  pm  wavelength  appears  to  be  more  effective 
than  the  A=4.875  pm  wavelength.  In  addition,  the  lower  PRF 
(18.7  MHz)  seems  to  contribute  to  higher  penetration  rates 
due  to  the  higher  energy  per  pulse.  It  is  recommended  in 
future  experiments  to  keep  the  pulse  energy  constant,  but 
increase  the  repetition  rate,  to  see  how  the  burn-through 
rate  is  affected.  It  would  be  useful  also  to  vary  the 
wavelength,  while  keeping  the  PRF  and  the  energy  per  pulse 
constant . 
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